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A class of integrable 2-dim classical systems with integrals of motion of fourth 

Qs, [ order in momenta is obtained from the quantum analogues with the help of 

ON ■ 

• deformed SUSY algebra. With similar technique a new class of potentials 



connected with Lax method is found which provides the integrability of 
corresponding 2-dim hamiltonian systems. In addition, some integrable 2- 



Oh 

<: 

' dim systems with potentials expressed in elliptic functions are explored. 

(N : 1. Introduction 

> 

o 

§ I Construction of classical integrable systems with additional integrals of motion is of consid- 
erable interest in Mathematical Physics (see |I| and references therein). Multidimensional 
integrable systems play an important role to describe the dynamics analogously to 1-dim 
manifestly integrable systems. In particular, they may serve as zero-approximations of per- 
turbation theory in the case of weak, nonintegrable perturbations. Variety of traditional 

O 

approaches to this problem exists starting from J.Kepler, S.Kowalewski till the Lax method. 

On the other hand a modern viewpoint on how to build classical integrable systems is 

based on the symmetries of related quantum systems 0. Recently the method for search- 
ed ' 

ing quantum integrable 2-dim systems was developed ||,[U with the help of a deformed 
supersymmetry (SUSY) algebra formed by intertwining differential operators of finite order. 
Supersymmetry 0-0, i.e. the construction of the isospectral pair of Hamiltonians, was 
proved [Q] to be in one-to-one correspondence to integrability of both Hamiltonians, i.e. to 
existence of a differential symmetry operator, which is polynomial in derivatives and which 
transforms solutions of 2-dim Schrodinger equation into other solutions with the same energy. 
Quasiclassical reduction of the deformed SUSY algebra [§] gave the factorization of classical 
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integrals of motion for the corresponding Hamiltonians ||. As a result, the structure of 
analytically resolved integrals of motion became more clear, and new classes of integrable 
potentials were found |§,||. 

In our paper we continue our study |9j of classical systems which integrability is induced 
or inspired by a deformed SUSY algebra for the relevant quantum systems. The consise basic 
construction of systems possessing a dynamical symmetry with the help of Higher derivative 
SUSY algebra is essentially supplemented with algorithms of searching analytical solutions of 
related nonlinear equations for coefficients functions of symmetry operators and potentials. 

In Sect. 2 the intertwining relations between a pair of quantum Schrodinger-type Hamil- 
tonians by general differential operators q^ of second order are investigated. The class of 
particular solutions of these relations is constructed for the cases of hyperbolic (Lorentz) 
gik = diag(l, — 1) and degenerate = (1,0) metric structures of operators q ± in second 
derivatives. The differential operators of fourth order in derivatives, which are symmetry 
operators for intertwined Hamiltonians, are built. In Sect. 3 the classical limit H — > for the 
Hamiltonians is considered, and the class of systems with integrals of motion of fourth order 
in momenta is obtained. In Section 4 a new class of integrable systems with potentials con- 
nected to the Lax method is derived using ansatzes and technique taken from 2-dim SUSY 
Quantum Mechanics (Section 2). Last Section 5 is devoted to description of some integrable 
systems expressed in elliptic functions. We stress that quite a few of the obtained potentials 
do not allow the separation of variables in known coordinate systems and some of them sofar 
have not been found. 



2. Quantum integrable 2-dim systems. 

In the 2-dimensional generalization ||, ||, f|] of Higher Order SUSY Quantum Mechanics 



TOU the intertwining relations of second order in derivatives are most essential: 



#<V = ? + # (2) , q-HM = H®q-, (1) 
H® = -h 2 A + V(x), A = 01 + 31 di = d/dxi, 
Q + = (<r) f = h 2 g ik (x)did k + hCi(x, h)di + B(x, h), 



where fi is the Plank constant and all coefficient functions are real. 
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It means that, up to zero modes of q ± , spectra of coincide and their eigenfunctions: 

^( 2 ) „ ~ g + ^ (2) . (2) 

The intertwining relations ([!]) lead to existence of the symmetry operators for 
the Hamiltonians H^ 2 \ correspondingly: 

#«] = 0, = q + q-, R {2) = q-q + , 1 = 1,2. (3) 

In 1-dim case [L(| analogous symmetry operators Ry> become polynomials of with 
constant coefficients. The distinguishing peculiarity of 2-dim case is existence Q of nontrivial 
dynamical symmetry operators R^ which are not reduced to functions of the Hamiltonians 
H®. 

It was shown in || that for the unit metrics = 5i k operators R^ can be written as 
second order differential operators (up to a function of ifw) and corresponding quantum 
systems allow separation of variables in parabolic, elliptic or polar coordinate systems. For 
all other metrics operators i?W are of fourth order in derivatives. 

The intertwining relations (P are equivalent to the following system of differential equa- 
tions: 

hd t C k + hd k Ci + h 2 Ag lk - (V™ - V^)g lk = 0; 

ft 2 AC* + 2hd t B + 2hg tk d k VW - (yd) - V^Q = 0; (4) 
h 2 AB + h 2 g tk d k d t V {2) + hCAV® - (V w - V {2) )B = 0. 

where the metrics go- is a quadratic polynomial in x%, x 2 . 

g u = ax\ + a x x 2 + h] g 2 2 = ax\ + a 2 x 1 + b 2 ; gvz = -^(2ax\x 2 + a\X\ + a 2 x 2 ) + b 3 . 

I. For the supercharges with Lorentz metrics (gi k = diag(l, —1)) : 

q+ =h 2 (dl-d 2 )+hC k d k + B, (5) 

a solution of can be reduced j|] to a solution of the system: 

5_(C_F) = -d + (C + F); (6) 
d 2 + F = d 2 _F, (7) 



where C\ =F C 2 = C±(x±) depend only on x±, respectively. Eq.(|7|) means that the function 
F can be represented as a sum F = F^(x + + xJ) + F 2 (x + — xJ). The potentials V^ l,2 > and 
the function B are expressed in terms of solutions of system (|]), ([F|): 

^ = ±^(C' + + C'_) + l(C 2 + + C 2 _) + ±(F 2 (x + -xj-F 1 (x + + x„)), (8) 

B = i(C+C_ + F 1 (x + + x„) + F 2 (x + -x_)). (9) 

The solutions for functions F, which admit additionally the factorization F = F + (x + ) ■ 
F-(x-), were found in In the present paper other solutions of (||) - ([7]) will be built. 



1) After substitution of the general solution of 



L[ 1^- /^)/(C + C_), (!()) 



C + J C- 

into (0), we obtain the functional-differential equation for functions L and A' ± = 1/C±(x±) : 
/A"' A"'\ 

- jr)HA + - A.) + 3« + - A.) + (X? - A' 2 )L"(A + - A.) = 0, (11) 

where V denotes the derivative of L with respect to its argument. Eq.(|TT|) can be easily 
solved for functions A± such that A'± = \ 2 A±, Then 

L{A + - A_) = a{A + - A_)- 2 + (3, 

for A± = a±exp(\x±) + 5±exp(—\x±) with a + ■ 5 + = ct_ ■ o~_ and a, (3— real constants. For 
A 2 > we obtain (up to an arbitrary shift in x±) two solutions: 

la) A± = k sinh(Aa;±), 16) A± = k cosh(As±). 

Then (|10|) leads to: 

la) F t (2s) = F 2 (2x) = A- - + k 2 cosh(2Ax), C± = k ^ 0, (12) 

cosh (Ax) cosh(Ax ± ) 

16) Fi(2s) = -F 2 (2x) = , A- . +fc 2 sinh 2 (Ax), C±= . - f. fc ^ 0. (13) 

smh (Ax) smh(Ax±) 

For A 2 < hyperbolic functions must be substituted by trigonometric ones. 
At last, in the limiting case of A = the solutions have the form: 

FA2x) = -F 2 (2x) =k lX - 2 + k 2 x 2 , C± = — , k ^ 0, (14) 

x± 

k 

Fi(2x) = -F 2 (2x) = hx 2 + k 2 x\ C± = ± — , k ^ 0. (15) 

x ± 



2) To find another class of solutions of the system (H), (0) it is useful to replace in (|i~0D C± 
by f±, such that C± = ±f±/ f' ± . Then F in ( |i~0D is represented in the form F = U (f+f-)f+fL 



with an arbitrary function {7. After substitution in ([/]) one obtains the equation: 

U+f- - flf' 2 )u"(f) + 3/0± - ^)c/'(/) + (jr - j-) W) = o, / = /+/- 

One can check that f± = a±exp(Xx±) + j3±exp(—\x±) and U — a + 46/ + /_ are its 
particular solutions (a, 6— real constants). Then functions 

= exp(Ax) + exp(— Xx)) + k 2 (a 2 + a 2 _ exp(2Ax) + (3\0t exp(— 2Xx j), 

—F 2 (x) = ki(a + P- exp(Ax) + (3+a- exp(— Ax)) + /c 2 (a+/? 2 exp(2Ax) + /3 2 a 2 exp(— 2Ax)), 
c ^ _ ± a± exp(Ax ± ) + j3± exp(-Ax ± ) 
\(ot± exp(Ax±) — (3± exp(— Xx±)) 

are real solutions of the system (j^j, (|7|) if a±, f3± are real for the case A 2 > and a± = (3± 
for the case A 2 < 0. 

3) To find a third class of solutions it is useful to rewrite @ in terms of the variables 

2(F 1 (x 1 ) + F 2 (x 2 ))d 1 (C + + C_) + Fi( Xl )(C + + C_) + F^(x 2 )(C + - C_) = 0. 

Its solutions are: 

3a) C + (x + ) = o\o 2 exp(Ax+) + 5i5 2 exp(— Xx + ) + c, 
C_(x_) = (TiS 2 exp(Ax_) + 020i exp(— Ax_) + c, 

F 1 (x 1 ) = 0, F 2 (x 2 ) = - 1 ; (17) 

(cr 2 exp(Ax 2 ) - o 2 exp(-Ax 2 )) J 

46 2 

36) C+(x) = C_(z) = ax 2 + c, Fi(xi) = 0, F 2 {x 2 ) = — ; (18) 

x 2 

3c) C + (:r + ) = a"ier 2 exp(Ax + ) + 5i5 2 exp(— Xx + ), 
C_(x_) = (7i5 2 exp(Ax_) + a 2 o"i exp(— XxJ), 

Fi, 2 (xi j2 ) = — ^| — — ± 7. (19) 

(cri )2 exp(Axi >2 ) ± 01,2 exp(-Axi >2 )) 2 

Let us remark that two additional solutions, analogous to 3a) and 3b), can be obtained by 
replacing Fi(x\) with F 2 (x 2 ) and vice versa. 
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After inserting these solutions (|12|) - fll9|) into the general formulae for potentials we 
obtain, correspondingly, the following expressions for potentials (^0|)-(^7|): 



^(1,2) 



T- 

1 
4 



sinh(Ax_| 



cosh 2 (Ax H 



+ 



sinh(Ax_ 



cosh 2 (Ax 



+ 



+ 



1 



cosh (AX2) cosh (Axi) 



cosh 2 (Ax + ) cosh 2 (Ax- 
k 2 cosh(2Ax2) — k 2 cosh(2Axi) 



+ 



(20) 



cosh(Ax + ) cosh(Ax_) 
sinh 2 (Ax + ) sinh 2 (Ax- 



k 2 

+ TT 



1 



sinh 2 (Ax + ) sinh 2 (Ax- 



1 



ki k\ 
. - 2/ . — r + . : 2 ._ — r + k 2 cosh(2Axi) + k 2 cosh(2Ax 2 ) 
.sinn {ax 2 J smh (Axi) 



(21) 



hk ( 1 



T 2 [xl + x 2 _ 



k' 2 



2 1 
X X 



—5 + -j + ^2(^1 + ^2) 



Let us note that the potential fl2~l|) with k 2 = and the potential 
the literature (c.f. for example 



(22) 

were investigated in 



=F 



tikf 1 1 \ A; 2 / 1 1 \ 1 



2 Vx5 



xl/ 8 Vxfj_ xl/ 4 
2a+/3 + (l =f 8/iA 2 ) + a 2 exp(2Ax+) + (3% exp(-2Ax 



fci(x 2 + x 2 ) + k 2 {x\ + x 2 ) 
+ 



y(l,2) 

8A 2 (a + exp(Ax + ) — (3 + exp(— Ax + )) 2 
2a_/3_(l±8frA 2 ) +a 2 _ exp(2Ax_) + (3 2 _ exp(-2Ax_) 
8A 2 (a_ exp(Ax_) — /?_ exp(— Ax_)) 2 

fci(a + /3_ exp(2Ax2) + exp(— 2AX2)) + k 2 {a 2 + (5 2 _ exp(4Ax2) + 

a 2 /3+ exp(— 4Ax 2 )) + k\{a + a- exp(2Axi) + (3 + (3- exp(— 2Axi)) + 

k 2 (a 2 + a 2 _ exp(4Axi) + (3 2 + (3 2 _ exp(-4Axi)) ; 

y(i,2) _ ±^( ai exp(Axi) - 81 exp(-Axi))(cr 2 exp(Ax 2 ) + 5 2 exp(-Ax 2 )) + 
(a 2 exp(Axi) + 5 2 exp(— Axi))(o"2 exp(Ax2) + 5 2 exp(— AX2)) + 



(23) 



(24) 



2c(o"i exp(Axi) — <5i exp(— Axi))(cr 2 exp(Ax 2 ) — <5 2 exp(— Ax 2 )) 
1 



+ 



4(<t 2 exp(2Ax 2 ) - 5 2 exp(-2Ax 2 ))' 



1 b 2 
y(i,2) _ ±2h aXl + -\a 2 (x{ + x\ + 6x 2 x 2 ) + ac{x\ + x 2 2 )\ + — ; 

4 x 2 

foA 

y(i>2) = ±— (cti exp(Axi) - 5i exp(-Axi))(cr 2 exp(Ax 2 ) + <5 2 exp(-Ax 2 )) + 



(25) 
(26) 
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g(°"i exp(Axi) + 8\ exp(— Axi))(af exp(Ax2) + 8 2 exp(— XX2)) + 

(27) 



{a 2 exp(2Ax2) — 82 exp(— 2AX2)) 2 (<7i exp(2Axi) + 5\ exp(— 2Axi)) 2 

The Hamiltonians with potentials (|20| ) - (|27j ) possess the symmetry operators 
q + q~,R( 2 ^ = q~q + , where can be obtained inserting solutions fljjQ - (|T9|) into fl5|). 



It is necessary to note that there are some singular points of potentials fl2~ID-(27) on 
the plane (xi,^)- Therefore both the asymptotics of corresponding wave functions and 
their behaviour under the action of supertransformation operators q^ (|5|) and of symmetry 
operators have to be investigated. In particular, for the potentials (|23|) - (pTD 

the operators q ± the preserve asymptotics of wave functions in the singular points, so the 
operators Rp-> , R^ are physical symmetry operators for these systems. For the potentials 
([HI), symmetry properties were discussed before (see for example, [|Tll ). 
II. For the supercharges with degenerate metrics = diag(l, 0) : 

q + = h 2 dl + hC k d k + B (28) 

Eqs.flf) lead to: 

C x {x) = -x 2 F[(xi) + Gi(xi); C 2 (x) = Fi(xi); 

V {1) = h(2G[ - x 2 F?) + ^x 2 2 (F 2 )" - x 2 (F 1 G 1 )' + Kfa) + K 2 (x 2 ); 

V {2) = hx 2 F'l + ^x 2 2 (F 2 )" - x 2 (F 1 G 1 )' + + K 2 (x 2 ); 

B = --{&! + x 2 F[') + l -G\ - \x\F x F>; + x 2 F 1 G[ - K^x,), (29) 

where the real functions Fi(xi), Gi(x\), Ki{xi) are solutions of the following system: 
h 2 C" ti 

1 + ^((Gl)" + 2Gf ) + G 1 K[ + 2G' X K X - F X {F X G X )' - G\G\ = m^; (30) 



h 2 



2 2 



—F± 1V> - h(F[G'{ + 2G[F[') - G 1 (2G' 1 F{ + G'[F, 



F[K[ + ^(F 2 )" - 2G' 2 F 1 - 2F{'K 1 = m 2 F 1 ; (31) 
^Gi(F 2 )"' + F{(F 1 G 1 )" + ZG\F X F« = m 3 F 1 ; (32) 
^F[(F 2 y" + F 1 F[' 2 = m 4 F 1 , (33) 
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and K 2 (x 2 ) is the polynomial of x 2 with constant coefficients: 

K 2 (x 2 ) = m - m x x 2 - -m 2 x\ - ^m 3 x 3 2 + ^m 4 x 2 . 

Several particular solutions of Eg . (|33D can be found. The constant function F\ = k\ is 
the solution of (|33|) for m 4 = 0. To find other solutions we define the new function U(F\) : 

U(F 1 ) = F[( Xl ), (34) 

to decrease the order of the differential eq. (fffi): 

V + ^U- + i-[7< - ^ = 0. (35) 



Inserting its known solution |T7J into the following equation for Fi(xi) is obtained: 



F x 4 + nF 2 + k)~ 1/4 dF 1 = Xl , n,k = Const. (36) 



m 4 

The integral fl3"6"|) can be written as a finite combination of elementary functions only 
in the case when two of three constants m 4 ,n, k are zero. Thus the solutions of Eq.(^) in 
elementary functions are: F\ = ki;F\ = x\jn\F\ = (3/2) 2 / 3 k 1 / 6 x^ 3 ; F\ = m 1 J 2 x\/4:. Below, 
for simplicity, we shall consider the solutions with particular values of constants 7714,71, A;, 
while solutions with arbitrary values of these constants will differ by some of coefficients 
only. To solve Eqs.(^) - ( p2[) it is useful to consider separately two cases: G\ = and 
G\ ^ 0. In both cases solutions with F\ = x± lead to potentials V^ 1 ' 2 ^ with separation of 
variables. Below on such solutions will be ignored. 

1) G x = 0. 

In this case the potentials for F\ = k\ correspond again to Schrodinger equations 

with separation of variables. More interesting choices F\ = x\ and F\ = x^ 3 lead, respec- 
tively, to potentials (I— arbitrary real constant): 

= T2hx 2 + Ix^ 2 + h^x\ + 6xlx 2 2 + 8xi) - ^{x 2 + 4x 2 2 ); (37) 
V™ = 7 Jl x ~ 2 ± | X2X -/3 + l,- 2 /3 (x 2 + 9 x?)9/x?/3 _ rr^ {9x2 + 4xi) _ (3g) 

2) G, ^ 0. 

2a) For F\ = k\ 7^ Eq . fl3T|) leads to the following equation for G\{x\) : 

G\dGi 



k — \m 2 G\ 



x x , 



which has the solutions in terms of elementary functions in two cases: when k > 0, m 2 
or k = 0,m 2 < 0. For the first one, after redefinition of constants and translation in x 2 , 



5h 9 hko , \ - 

X? + -^(l±l) Xl 



36 

2 3/cimi 



2/3 k\k 2 _2/3 . 

—x 2 x x + 



. 2/3 
)X{ 



k\ 

mix 2 + —■ 



(39) 



The second case (k = 0, m 2 < 0) leads to separation of variables. 



2b) If F x 



x 



2/3 



and F\ = x\ the general solutions of Eq.(^) can be found and after 



substitution into (|30|) , (|3l| ) give the function K\(x\). Corresponding potentials are: 



^(1,2) 



7/T _ 2 2ft 
36 1 9 



3mi 

8*r 



,2/3 



+ 



4 



XT 



- m x x 2 + — (1 + 



2/3 (^ + 



5fcl 2/3 . 
— ^2^1 + 



15/ci)x: 



2- 



^(1,2) 



3ft 2 k 2 1 

( — =F %ki + y)^r 2 =F 2ftx 2 + ~(x\ + 6x 2 x 2 + 8x 2 ) 



-(12k 2 + m 2 )(x\ + 4x 2 



fcm 2 + 6fc 3 + 4/ci 



-x 2 . 



(40) 



(41) 



Because all potentials (|3~?D - ( STJ ) are singular at the x 2 axis, similar to the case of Lorentz 
metrics, it is necessary to investigate separately the behaviour of wave functions at x\ — > 
0. Straightforward though cumbersome calculations show that the fourth order differential 
operators (see (Bp) preserve the asymptotics of wave functions for systems fl3"7|) - 

((41|) and play the role of true dynamical symmetry operators. 



3. Construction of 2-dim integrable classical systems 
by the limit h — ► 

Quantum dynamical symmetries which were found by the intertwining method in Sect. 2 have 
their natural analogs - integrals of motion - in the corresponding classical systems. These 
integrals of motion are polynomials of fourth order in momenta. Similar to Sect. 2, it is useful 
to consider separately the classical limits for Lorentz and degenerate metrics. 

I. For the Lorentz metrics in the limit fi — > classical supercharge functions become: 

<ld = ~4p+p- ±i(C_(x_)p+ + C + (x + )p_) + 5(x_,x+). 
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From Eqs.(|8D and ([3]) we find that the classical Hamiltonian 

h d = 2(p 2 + + pi) + i(Cj + C 2 _ ) + \[F 2 {x + - x_) - Fi(x+ + (42) 
has the additional integral of motion: 

I = I6p 2 + p 2 _ + C\p 2 _ + C 2 _p\ - 2{F X + F 2 )p + p_ + B 2 . (43) 



Such sort of classical systems was considered in the literature (see [T^] and references 
therein). Usually the functional equation, which provides existence of integrals of motion for 
corresponding classical systems, is solved by the Lax method. Comparison of (42) and ( 43|) 



with notations in [13] leads to the relations: 



v x = l/8C 2 Jx 1 ),v 2 = l/8C 2 _(x 2 ),v 3 (x_) = l/16F 2 (x^,v 4 (x + ) = -l/16F 1 (x + ), (44) 



and functions Vk must satisfy [13|] a certain functional equation (see Eq.(f4"T|) below). In the 
next Section we prove that (f44|), where C±,Fi )2 are solutions of the system @, ([?]), satisfy 
also Eq . (|47|) . Moreover some additional solutions of this equation will be found in Sect. 4. 

II. Let us study the integrable classical systems which can be obtained in the limit 
h — > from SSQM systems with degenerate metrics in q^. The classical supercharges have 
the form q d = (q~i)* = —p\ + iCk(x)pk + B(x) and the Hamiltonian 

Hd =pl + \x 2 2 d 2 x Fl - x 2 {F x G x )' + K x {x x ) + K 2 (x 2 ) (45) 

has the integral of motion of fourth order in momenta: 

/ = QtiQd =Pt+ (C? " 2B)pI + C 2 p 2 + 2C x C 2 p x p 2 + B 2 . (46) 

All functions in (^5]), (^) were defined in the previous Section, where we have to put h = 0. 

Thus in the case of degenerate metrics the following classical integrable systems are 
obtained (new definitions of constants were used for some of these systems): 

1) V = ^-(x{ + <ox\x\ + 8x2) + m{x\ + 4^2) + lx x 2 , 

I = p x + (6x^2 + mx\ + x\ + 2lx x 2 )p\ + x\p\ — 

x 4 

Ax\x 2 p x p 2 + {x\x\ + mx\ + + ^r 2 ) 2 ! 
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2) V = \ Xl 2, \\x{ + x 2 ) + m(9xl + 4x 2 ) + kx 2 / 3 , 

I — Pi + {-x\xi 2 ^ + 2rex 2 ^ 3 + 18mx 2 + Xi^)p\ + x^ 3 p\ — 
9 

-x 2 x\ /3 pip 2 + i^xlxi 2/3 - rex 2 / 3 - 9mx 2 - ^Xi /3 ) 2 ; 

3) V = -^^ L X 2 Xi 2lZ + \{k\ + ^ mi )x 2 1 /3 - 7711X2, 

3 4 re 2 

r 4 , An2 i 3fci7?7i 2 /3 2^!^ _ 2 /3 . , 2i 2 , 1 2 2 , 

1 = Vi + [^{k 2 + fc )Xi ^~x 2 x 1 +k 1 \p 1 + k 1 p 2 + 

2k 1 k 2 x{ /3 p 1 p 2 + [^^x 2 x^ 2/3 + hjzl - 3 ^ 1 )xi /3 - y] 2 ; 

A\ T/ 1 - 2 / 3 r 9 2 , 2n 5fc 2/3 . 3m 2/3.^2 . k /, . , r , v 2 

4) V = -xi (-x x + x 2 ) - yx 2 xi 7 + gT-^i + -j^i - mx 2 + -(1 + 15re)x 2 , 

T 4 , 2 -2/3 10/C 2/3 k 2 2 , 2/3 . 4/3n 2 , 4/3 2, 

J = Pi + (g^i -— ^a;/ + y Xi + ^A rXi +X/ k + v p 2 + 

2/3/ 2 —1/3 \ / 1 2 —2/3 2/3 3?7ll 2 /3 1 4/3 ^ 2 2\2 

2x/ (— -x 2 Xi + kxi)pip 2 + (-x 2 Xi + rex 2 Xi' - -^-x/ - -x/ + — arj ; 
3 9 ore 2 4 

1 re 2 3 

5) V = ~(x? + Qxfxj + 8xt) + m{x\ + Ax 2 2 ) + -^x^ 2 - (h - 2mk - -re 3 )x 2 , 

k 2 

I = p{ + (6x 2 x 2 + 2(m + re 2 )x 2 - 2reix 2 - 4rex 2 x 2 + x{ + y^ 2 + rerei)jj 2 + 

re 2 1 k 2 kki 
x\p\ + 2xi(rei - 2x 2 1 x 2 )p 1 p 2 + {x\x\ + (m + y)^ 2 , + + ^ lX2 ~~ ^ x ^ 2 ' 



These potentials are not new: they were found by other methods in [I4"| - [IB|. 

In conclusion of this Section we formulate the procedure of construction of integrals of 
motion in terms of Classical Mechanics objects. For any classical Hamiltonians H c i and a 
complex function q^(x, p) = (<?^)*, polynomial in momenta, such that: 



{Qd, H d} = if{x,p)qti, {{qti)*,H d } = -if(x,p)(q, 



cl) ) 



with arbitrary real function f(x,p), the classical factorizable integral of motion I = q^ ■ q cl 
exists ({, }— Poisson brackets). 
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4. Integrable systems connected with the Lax method. 

Let us consider classical systems with potentials of the form: 

V(Xi,X 2 ) = fl(^l) + v 2 (x 2 ) + v 3 (xi - x 2 ) + v 4 (xi + x 2 ). 



It is known [I3~ that these systems have the integrals of motion of fourth order in momenta: 

1 = T^PiPl + v 2 (x 2 )pl - (v 3 (xi - x 2 ) - v 4 {xi + x 2 ))pip 2 + Ui(xi)p2 + /, 

if the functions V\,v 2 , v 3 , V4 satisfy the basic functional equation: 

[v^Xt + x 2 ) - v 3 (xi - x 2 )][v' 2 \x 2 ) - v'Kxx)} + 2[v'l(xi + x 2 ) - i#(a?i - x 2 )} ■ 
[v 2 (x 2 ) - Vi(xi)] + 3w4(xi + x 2 )[v' 2 (x 2 ) - v[(xi)] + 

3v' 3 (xi - x 2 )[v' 2 {x 2 ) + ui(xi)] = 0. (47) 

There is a list of known particular solutions of Eq . (|47|) in the book |13[ . To search for 
new solutions of this equation it is useful to rewrite it in the equivalent form: 



d 2 (vv 2 + 2v 2 d 2 v) = di(vv[ + 2vidiv), 



(4? 



where 

v = v±(xi + x 2 ) - v 3 (xi - x 2 ). 
The general solution of Eq.(f48l) is: 



(49) 



G 



dx\ 



dx-. 



+ L 



dx\ 



dx-, 



^i^O; v 2 ^0, (50) 



where G and L are arbitrary functions of their arguments. Thus the problem is reduced to 



searching for the functions v of the form (|50|) , which satisfy the condition 



{d\ - d 2 2 )v = 0. 



(51) 



In particular, it is easy to check that all solutions which were found from SSQM in K — > 



limit in the Sect. 3 (see Eq.(44)) are particular solutions of Eq.(^) and have the form ( |50| ) 
with G = 0. 

Let us apply the technique, which was used in investigation of the system fl6|), (0) in the 
framework of SSQM, to find the new particular solutions of eq. (4~8). 
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1) If the function v is factorizable v = Ui(xi) ■ 1x2(^2), Eq. 
variables and its solutions have the form: 

n k [a k exp(y/\x k ) + b k exp(-\/~\x k )] + l k 



V3 



(a fc exp(v / Ax fc ) - b k exp(-v / Ax fe )) 2 
a\b 2 exp(VA • xJ) + a 2 b\ exp(— VA • x_ 

aia 2 exp(V'X ■ x + ) + b\b 2 exp(— v^A • x A 



admits separation of 



k = 1,2. 



(52) 



where for A > all constants are real and for A < a k = b k and n k , l k are real. 
2) Let us introduce new functions 



V 2 



(53) 



Then Eq. 



|) means that (derivatives of G and L are taken in their arguments): 



( G + L \wi Wi 

{W[ 2 -W?){G" + L" 



+ 3 



(W» - w%)G' + {W'{ + W%)L' 



+ 



0. 



(54) 



It is possible to construct several particular solutions of (|5 
2a) 

Wi.,2 = o"i,2exp(Axi j2 ) + 5i,2 exp(— Axi ]2 ) 



(55) 



where constant A 2 is real and a k , S k are complex. If these constants satisfy the condition 



ai5\ = a 2 5 2 , 



(56) 



then substitution of (pq) into (p4|) leads to an equation with separable variables. Its solutions 
are: 



G(Wi + W 2 ) 
L(W 1 - W 2 ) 



(Wi + W 2 ) 2 
a 2 



+ a{Wi + W 2 f + fin 
+ a(W 1 - W 2 f + fa, 



(57) 
(58) 



with arbitrary constants a,cti,(3i(i = 1,2). 

Let us consider the case with A 2 > 0. Constants 0i,5i must be both real or both pos- 
itive because of requirement that functions v n (x n ) for n = 1,2 should be real. Analogous 
arguments work for the pair a 2 ,S 2 . The condition (|56|) for real (cij, 5a leads to two options: 



Wi(x) = W 2 (x) = fccosh(Ax), 

Wi(x) = W 2 (x) = fcsinh(Ax), k e R, 



(59) 
(60) 
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and the function v is real for real constants a, aj, /3i(i = 1, 2). 

In the case of real (a\, and imaginary (a 2 , 8 2 ) the solutions take the form: 



Wi = fcsinh(Axi), W 2 = ik cosh(Ax 2 ), (61) 
Wi = fccosh(Axi), W 2 = iksinh{Xx 2 ), k G R. (62) 

Let us remark that in this case the arguments of the functions G and L are complex conju- 
gated and in order to have a real function v the following relation must be fulfilled: 

L*(Wi + W 2 ) = -G(Wi + W 2 ). 



Therefore a 2 = -a\- (3 2 = a = a* in fl57|), (H). 

2b) 

Wx = W 2 = gx 2 ; g 2 G R. (63) 

It follows from Eq.(||) that: 

G = a + a 1 (W 1 + W 2 ) + a(W 1 + W 2 ) 2 ; L = b - - W 2 ) 2 + 



4 V " (Wi + iy 2 ) 2 ' 

where all constants are real. 



Thus the functions v n (n = 1, .., 4) for solutions (^) - (|63|) are, correspondingly: 

k 



vi(x) = v 2 (x) 



sinh 2 (Ax 



k\ k 

v 3 (x) = v 4 (x) = 2 / \ \ + . , 2/Xx , + h cosh(2Ax) + k 4 cosh(Ax); (64) 
smh (Ax) smh (^-) 



Vx{x) = v 2 (x) 



cosh (Ax 



k\ k 2 

v z( x ) = . \ + ■ , 2,\x\ + k 3 cosh(2Ax) + k 4 cosh(Ax) 

smh (Ax) smh (^) 



v *( x ) = ^Xf^ ~ ■ ^2,\ X \ + fc 3 cosh(2Ax) - k 4 cosh(Ax); (65) 
smh (Ax) smh {*§■) 



v i( x ) — j} — T) v 2 (x) = 27 — -; 

cosh (Ax) sinh (Ax) 

, . , . ki + k 2 sinh(Ax) . . , . , s , N 

y 3 (x) = w 4 (x) = -n— — h K3Cosh(2Ax) — /c4smh(Ax); (66) 

cosh (Ax) 
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vi(x) 



vAx 



vAx 



k ( 
sinh (Ax) 
k\ + k 2 sinh(Ax) 

cosh 2 (Ax) 
k\ — k 2 sinh(Ax) 

cosh 2 (Ax) 



k 



cosh 2 (Ax) ' 
+ k 3 cosh(2Ax) + ki sinh(Ax) 

+ k 3 cosh(2Ax) — ^4 sinh(Ax) 



vi(x) = v 2 (x) 



- kx 
9i 



-2 



VAX) 



VAx) 



V2 



k\X 2 + k 2 x 2 + k 3 x 4 + /C4X 6 , 

92 



(67) 

(68) 
(69) 



(Si exp(Ax) — 0i exp(— Ax)) 2 ' (52 exp(Ax) — a 2 exp(— Ax)) 2 ' 

V3(x) = ki(5i<7 2 exp(Ax) + ai5 2 exp(— Ax)) + k 2 (S 2 al exp(2Ax) + a\5 2 exp(— 2Ax)); 

vA[x) = ki(6i5 2 exp(Ax) + o\o 2 exp(— Ax)) + k 2 [a\a\ exp(2Ax) + 5 2 5\ exp(— 2Ax)). 

The solutions fl6"5|) - ([w]) are absent in the list of |L3| and to our knowledge they are 
novel. As to expressions (|69"D , they are present in |L3| only for a\5i > 0, 02^2 > 0. 

In conclusion, let us note that one can easily check some invariance properties of Eq . fl47|) . 
In particular, from arbitrary solutions (|52"D, (Q) - ( |5"9"| ) one derives new ones with: 

v 4 (2x) — > Vi(x); v 3 (2x) — > v 2 (x); vx(x) — > -v 3 (x); v 2 (x) — > -v 4 (x). (70) 

Eq.(^7|) is invariant if V\ >2 — > f i )2 + c and f 3j4 — > f 3i4 + c with arbitrary constants c, c. It is 
also invariant under dilatation of all arguments Xj — + Axj. 



5. Integrable systems with potentials, expressed in el- 
liptic functions. 

In this Section we formulate the method of construction of new integrable systems from 
solutions Vi of Eq. (pT7|) , which were found in previous Section. If we define new functions: 

Wf = /i(WU W? = f 2 (W 2 ), (71) 



Eq.(p4|) takes the form: 

[G(W 1 + W 2 ) - L(Wx - W 2 )][f 2 \W 2 ) - fliWx)] + 
2[G"(W 1 + W 2 ) - L"{Wx - W 2 )][f 2 (W 2 ) - hiWx)] + 

3G"(W! + W 2 )[f 2 {W 2 ) - f[{W x )\ + 3L'(Wi - W 2 )[f 2 {W 2 ) + f[(W x )\ = 0. (72) 
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This equation has the same structure as Eg. (|47|) , thus from ([71]) we obtain: 

W? = v 1 (W 1 )) W? = v 2 (W 2 ); 

G(W 1 + W 2 ) =v 4 (W 1 + W 2 )- L(W 1 -W 2 ) = -V3(Wi-W 2 ), (73) 

where we can use solutions (|52|) , (|64"D - (|69[) for v n (x) {n = 1, ...,4) with arguments 
Xi >2 replaced, respectively, by W\ j2 . The first pair of equations in (|73|) defines functions 
Wi(xi), W 2 (x 2 ) and last two equations define new functions G, L. Substituting these set of 
functions into (|53[) and (|50|) , we find some new solutions V n of the same Eq. (f47[) from already 
known solutions v n (see fl52|), (|64|) - (|6^)). 



Let us consider several examples of this method of reproducing new solutions. 



1) The first attempt to start our procedure from the simplest solutions (|68|) leads to 
discouraging result: we obtain the same solution. But we can firstly transform f 12 w 3 4 
in (|68|), using the invariance property ([T(]) mentioned at the very end of Sect. 4. As follows 
from (|73"D, the functions W± j2 (2:1,2) are defined from the equations (we omit indices i = 1,2): 

W >2 (x) = k x W~ 2 + k 2 W 2 + hW A + k A W 6 + k , 

with constant fcj. It is useful to rewrite them in terms of functions U(x) = ^W 2 (x) : 

U' 2 (x) = h + k U + k 2 U 2 + k 3 U 3 + hU\ (74) 

From (|50|) and ( ^3l) we obtain new solutions: 

y(xi) t/(x 2 ) 

Vl = ^Z7^)' ^ = ^f7^)' (75) 



y 4 ( Xl + X2 ) - V 3 ( Xl - X2 ) = g , TT ,\ v T ;,*L 2 , (76) 



^(^i)^^) 
(U( Xl ) - f/(x 2 )) 

where arbitrary constants g\,g appear because solutions v n are defined from Eq.(^7[) up to 
constant factors. We can check directly that for U(x), which satisfy Eq.([74]), r.h.s. of ( [76|) 
is the solution of Eq. (|5l|) . 

When r.h.s. polynomial in Eq. ([74]) has degenerate roots, U(x) can be expressed through 
elementary functions and corresponds to solutions fl52|), (|64D - (^), found above. When all 
roots are simple, U(x) can be given in terms of elliptic functions, 
la) U (x) = p(x) + b, 

where b is constant and p(x) is the Weierstrass function with semiperiods uii, u 2 (their values 
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depend on constants kj). Eqs.(^), (|?6| ) lead to new solutions: 

Vi(x) = V 2 (x) = a x p{x + uji) + a 2 p(x + uj 2 ) + a 3 p(x + (u l + lo 2 )/2); 
V 3 (x) = Vi{x) = ap(x), 

where a is arbitrary constant and constants ak(k = 1,2,3) satisfy the following condition: 

k ijtj 

lb) U (x) = snx+5, where b is constant and sn (x) is the Jacobi function with modulus 
k (it depends again on constants ki). In this case new solutions are: 

sni + b 



VAx) = V 2 (x) 



cn 2 x ■ dn x 



2 ' 



V 3 (x) = c( J j- - fc 2 cn 2 (x/2) - k 2 ), V 4 (x) = c(l - k 2 )(- ' ' 



"sn 2 (x/2) ' cn 2 (x/2) sn 2 (x/2) 

lc) [/ (x) = dnx + 6. 
Correspondingly, the new solution takes the form: 

Vi{x) = V 2 (x) = a J 2 ; V s (x) = V 4 (x) = c(^— - + (1 - A; 2 )cn 2 (x/2)). 
sn 2 x ■ dn x sm(x/2) 

2) The second solution of (f|7]), which we can take as the starting point of proposed 
procedure, is one of solutions in Eq . (|52]) : 



vi(x) = v 2 (x) = acosx + c, a > 0, c > 0, 

v 3 (x) = k 1 (sm(x/2)y 2 + £; 2 (sin(x/4))~ 2 , u 4 (x) = k 3 (sin(x/2)y 2 + £; 4 (sin(x/4))~ 2 . 
Then according to Eq. (|T3"D , functions W(x) must be found from equation: 

W' 2 (x) = a ■ cosW(x) + c, 

which solution can be expressed through the Jacobi function with modulus k {k 2 = 2a/(a + 

o)): i 

W(x) = arccos(l — 2(k ■ sny)~ 2 ), y = —yj (a + c) • x. (77) 



Thus the new solution of Eq. fl47D is: 

2, /«n cn 2 (y + /2) dn 2 (y + /2) dn 2 (y+/2) 

V 4 X! + x 2 = a lS n 2 y + /2 + a 2 2 ) y+, ) + a 3 ) + « 4 2 L • 

dn (2/4-/2) sn%+/2)cn% + /2) cn%+/2) 

2/ /.x cn 2 (W2) dn 2 (W2) 1 
V 3 (xi - x 2 ) = a 2 sn (j/_/2) + ai— — + a 3 ^- — — — — - + a 4 



V 1 {x) = V 2 (x) : , 
cn z y 



dnV/2) ^sn 2 (y_/2)cn 2 (y_/2) %n 2 (y_/2)' 

a 
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where do, CLk are arbitrary constants. 

In conclusion, we mention briefly the analogous method of construction of new integrable 
systems in Quantum Mechanics (see Sect. 2). In this case the procedure is based on Eq.(p|) 
of Sect. 2. Similarly to (|7T|), we introduce in ( |TTD new functions M±(A±) : 

A»(x+) = M + (A + ), A'_ 2 (x„)=M_(A_). (78) 

Then Eq. (pTTf) takes the form: 

[M»(A+) - M'L(AJ)]L{A + - A+) + 3[M' + (A + ) + M'_(A_)]L'(A + - A+) + 
2[M + {A + ) - M_(A_)]L"(A + - A+) = 0. 

The general solution of this equation can be found in |T7J and the corresponding functions 
M±(A±), L(A + — A J) can be used here to find A±(x) from (|78|). Substitution of these 
functions A±{x) into M±(A±), L(A + — A J) and Eq.fllOD leads to new solutions of the system 
Eqs.©, (0) for the quantum case. 
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